High-order Harmonic Generation in a Time-integrated Quantum Transition Picture 
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We develop a numerical scheme to investigate the high-order harmonic generation (HHG) in in- 
tense laser-matter interactions. Tracing the time evolution of every electronic laser-field-free state, 
we "observe" the HHG in a time-integrated quantum transition picture. Our full-quantum simu- 
lations reveal that continuum electrons with a broad energy distribution contribute equally to one 
harmonic and the excited state also plays an important role in the molecular HHG. These results 
imply a laser-intensity-dependent picture of intramolecular interference in the HHG. 

PACS numbers: 42.65.Ky, 32.80.Rm 
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Owing to the important applications in producing 
burst of high energy photon, high-order harmonic gen- 
eration (HHG) is a main focus of intense laser-matter 
physics 1]. In recent years, it is also shown that the 
HHG can be used to image molecular orbital0. The 
HHG can be well understood by a semiclassical recol- 
lision modelps']: (i) ionization of the active electron by 
tunnelling; (ii) acceleration of the electron in the laser 
field; (iii) recombination of the electron into the bound 
state to emit a high energy photon. Another widely used 
approach to describe the HHG is the strong field approxi- 
mation (SFA)[||, which can be regarded as the quantum- 
mechanism version of the semiclassical recollision model. 

Numerical investigation on the HHG of H^ demon- 
strates an interesting phenomenon: a minimum appears 
in the HHG spectrum that relates to the molecular ori- 
entation. Using a point-emitter model, the minimum 
is identified as arising due to intramolecular two-center 
interference and the position of the minimum is con- 
cluded not to rely on the laser intensity [H]- This min- 
imum has generat ed g reat theoretical and experimental 
interests [10, i, i, M,M, El El, [11 ■ Two experimental 
groups reported the observations of the minima in the 
HHG spectra of C02[3,[i|- They attributed the minima 
to two-center interference. However, the positions of the 
minima are different in their measurements (the 27th vs 
the 33rd orders). Other than interference, the theoretical 
calculations proposed that the ground state depletion is 
also a possible mechanism for the different modulations 
observed 9] . Zhou et al showed that two-center interfer- 
ence is responsible for the minima in the HHG spectra 
of CO2 measured in their experiments [T^ . Nevertheless, 
they found that the interference pattern is subjected to 
the Coulomb effects. Very recently, the experimental ob- 
servations of the interference minima in the HHG spec- 
tra of H2 were reported[14]. But the minima are found 
to be related to the laser intensity. This relation is at- 
tributed to the motion of the parent nuclei there. All 
these strongly call for a revisit of the HHG, especially, a 



quantitative description of intramolecular interference in 
the HHG. 

The core of the HHG is the recombination process and 
the recombination is an intermediate process that can 
not be observed experimentally. In addition, the semi- 
classical models that omit the Coulomb effects can not 
describe the HHG quantitatively. Thus, in this paper, 
we develop a numerical scheme to investigate the process. 
We directly solve the Born-Oppenheimer time-dependent 
Schrodinger equation (TDSE), and we project the wave- 
function on the eigen-states of the field-free Hamiltonian 
and trace its temporal evolution. With the method, we 
achieve a time-integrated quantum transition picture of 
the HHG. It gives us direct information about the elabo- 
rate role of every electronic state in the HHG and allows 
us to identify the whole energy transfer among harmon- 
ics, the laser field and the electrons in the HHG. Our 
simulations reveal that the major contributions to one 
harmonic come from continuum electrons with a broad 
energy distribution and the first excited state also has a 
significant contribution to the molecular HHG. As a re- 
sult, the position of the interference minimum is affected 
by the laser intensity in the frozen nuclei case. 

Below, our discussion is first made for the ID case, we 
then validate the main results to the 2D and 3D cases. 
The Hamiltonian of model molecules H^ or hydrogen- 
like atoms, all with the ionization potential Ip = 1.11 
a.u. studied here, is H{t) = -h V{r) - r • E(t) (in 
a.u. h ^ e = rrie = 1). Here, V{r) is the Coulomb poten- 
tial. In the ID and 2D cases, we use the soft-Coulomb 
potential V(x) = —. = \ ■ = and 

Vix, y) = —, -I- 

tively, where Z is the effective charge, R is the internu- 
clear separation. For H^, R = 2 a.u.. R = a.u. corre- 
sponds to the hydrogen-like atom. In the case of the 3D 
atom, it is V^(r) = —Z/r. E(t) — e£ sinLUgt is the exter- 
nal electric field with the amplitude £ and the frequency 
ujQ. e is the unit vector along the laser polarization. We 
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Figure 1: Harmonic spectra of molecules [(a),(c),(e),(f)] 
and hydrogen- like atoms [(b), (d)] calculated using different 
methods. The black curves denote the results obtained using 
Eq. f , the red and green curves denote those obtained using 
Eq. 2 with n' = 0, f and n' = Q, respectively. For the 2D and 
3D cases in (d)-(f), also see the context for the details. 

use 6 to denote the angle between the molecular axis and 
the laser polarization. Our calculation will be performed 
for 780 nm trapezoidally shaped laser pulses with a total 
duration of 10 optical cycles and linear ramps of three op- 
tical cycles. The details of solving the ID and 2D TDSE 
can be found in Ref.[il]. The 3D TDSE is solved follow- 
ing the program described in Ref . [lB| • The coherent part 
of the HHG spectrum is obtained byQ 

F(c.) = j (VW|g-VV|^(i))e'"*di, (1) 

where \'ip{t)) is the time-dependent wavefunction of H{t), 
and Lo is the emitted proton frequency. 

To achieve a full-quantum description of the recombi- 
nation in the ID case, we project the wavefunction {ipit)) 
on the eigen-states of the field-free Hamiltonian Hq = 
py2 + V{x), that is \ij{t)) = En<^n{t)\n)+J dpcp(t)|p). 
Here, ipdi^jt) — J c?pcp(t)|p) represents the continuum 
electronic wave packet. |n) is the bound state of Hq with 
negative eigen-energy En', |p) represents the continuum 
state of Hq that has positive eigen-energy denoted by 
Ep. The eigen-states and eigen- values are obtained by 
numerically diagonalizing the Hamiltonian Hq. 

The recombination in the HHG can be simulated by 
the evaluation of the dipole movement between the " con- 
cerned" bound states and the continuum electronic wave 
packet i/.c(r,t), i.e., D{t) = E„' <' Wi^'lV^IV^clr, t)). 
Here, n' indicates the concerned bound state that the 
continuum electrons can transit back to. Because the 



Figure 2: The contribution of |G„'(p,cj)p to certain single 
harmonic orders [(a)-(d)] and the related transition dipoles 
{n'\W\p) [(e)-(f)] for the spectra in Fig. 1(a) and (b). The 
red curves denote the results obtained with n' = 0, the black 
curves denote those obtained with n' = 1. 

population of other higher bound states is found to be 
negligible in our simulations. In this paper, we focus on 
the dipole movement in the cases of n' = and n' = 0, 1, 
i.e., the dipole movement relating to the transition back 
to the ground state |0) and to both the ground state |0) 
and the first excited state |1). The Coulomb effects, the 
ground state depletion, and the contribution of higher 
bound states those are omitted in the SFA, have been 
accurately incorporated into the expression of D{t). 

The Fourier transformation of D{t) gives the complex 
amplitude describing the HHG at the frequency u 

F{lj)= J dtD{t)e"^' ^ J dpG{p,Lu), (2) 

where G(p, uj) = J2n' (P' ^) 

Gn'{p,uj)^ {n'\WV\p) J dta:At)e''^'cp{t). (3) 

Here, (n'|VF|p) is the transition dipole matrix element 
between the continuum |p) and the bound state \n'). 

Eq. 2 shows that the contributions to one harmonic 
UJ come from the "returned" continuum electrons with 
varied energy Ep, each of which is weighted by the prob- 
ability amplitude G{p,u!). The HHG picture revealed 
here is somewhat different from the semiclassical ones. 
The electronic momentum in the semiclassical models is 
the instant momentum that is modulated by the laser 
field. The momentum in Eq. 2 is the " laser-field- free" 
momentum. Thus, Eq. 2 implies a "static" quantum 
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transition picture of the HHG, i.e., the transition of the 
continuum electron |p) with the time- integrated amph- 
tude G„'(p,Ci;) back to the bound state \n'). From the 
"static" picture, we can read some important messages 
of the recombination, as shown below. 

In Fig. 1(a) and (b), we compare the harmonic spec- 
tra of ID and its reference atom, calculated using 
Eq. 1 and Eq. 2, respectively. We expect that the com- 
parison can give an insight into the mechanism of the 
hhg[i1]. First, we see the black curve in Fig. 1(a) for 
shows a pronounced suppressed region, i.e., a broad 
hollow around the 27th order indicated by a dotted ar- 
row. This striking hollow structure is a two-center in- 
terference characteristic of the molecular HHG spectrum 
and will be discussed in detail below. In addition, the 
red curves in Fig. 1(a) and (b) basically reproduce the 
behavior of the corresponding black curves in the plateau 
region. This verifies the applicability of Eq. 2. The dif- 
ference in the low energy region is due to the omission 
in Eq. 2, the bound-bound transitions. Especially, in 
Fig. 1(a), the green curve obtained with n' — Q (i.e., 
ignore the first excited state) shows a deeper hollow than 
the red one that includes the first excited state. In Fig. 
1(b) of the atomic case, the two curves are also analo- 
gous. Thus, Fig. 1(a) and (b) reveal a main result of this 
present paper that the first excited state, which is usu- 
ally omitted in the HHG, could play an important role 
in the molecular HHG in the hollow regime. 

Next, we explore the quantum mechanism behind these 
phenomena. Fig. 2 plots the contribution of |G„'(p, w)p 
to certain single harmonics and the corresponding tran- 
sition dipoles (ri'|Vy|p) for the spectra in Fig. 1(a) and 
(b). We choose two typical harmonic orders, the 27th and 
69th, that locate at the interference regime and plateau 
regime of the HHG spectrum in Fig. 1(a), respectively. 
Fig. 2(a)-(d) obviously show that for the atomic cases, 
the amplitude is mainly contributed from the electrons 
that return to the ground state. For the molecular cases, 
electrons that transit back to either the ground state or 
the first excited state could contribute the amplitude sig- 
nificantly, especially in the low energy regime. 

The red curves of |G'o(p,ti')P in Fig. 2(a) and (c) of 
the molecular cases show a dip structure around energy 
E{tp) = 0.35 a.u., indicated by a blue arrow. The black 
curves of |Gi(p,a;)p show similar structures except that 
the position of the dip structure shifts to a higher elec- 
tronic energy of 2.25 a.u.. Because its amplitude is small, 
we can safely ignore its contribution in the following dis- 
cussion. The first dip that locates at -E(p) = 0.35 a.u. 
is important and corresponds to the pole of the molecu- 
lar transition dipole in Fig. 2(e). The latter is believed 
to be relevant to two-center interference of the diatomic 
molecules [l^ [l3|- As we integrate the amplitude over 
electron energy, this "two-center interference effect" re- 
mains visible for the case of the 27th, whereas it vanishes 
for the case of the 69th, as presented in Fig. 3(a)-(d). 

In the left panels of Fig. 3, we plot the integrated con- 
tribution of I !o d'pGn'(ji,i^)\'^. For the 27th order. 
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Figure 3: The contribution of | dpG(p,a')p to certain sin- 
gle harmonic orders for the spectra in Fig. l(a),(b) and (d). 
The red curves denote the results obtained with n' — 0, the 
black curves denote those obtained with n' = 0, 1. 



in Fig. 3(c) of the atomic case, the red curve shows that 
the primary contributions of | dpGo(p, w)p to this or- 
der come from those electrons with a broad energy distri- 
bution centered at i?p=0.467 a.u., as indicated by a blue 
arrow. This energy value agrees with the energy relation 
Lu = 27^0 = Ep + Ip {ujQ is the laser frequency) [1, [ll] . 
In addition, the black curve of rt' = 0, 1 almost super- 
poses the red one. In Fig. 3(a) of the molecular case, the 
red curve says that the primary contributions come from 
those electrons with energy smaller than i?p=0.467 a.u.. 
In Fig. 2(a) there is a dip in the red curve of |Go(p, uj)]"^ 



around E„ 



0.35 a.u.. We conclude that two-center 



interference responsible for the dip will induce a strong 
suppression on this order. Due to the suppression effect, 
the first excited state will also play an important role in 
this order. As we can see, the black curve is higher than 
the red curve here and this order in Fig. 1(a) is strongly 
suppressed. 

For the 69th order, both Fig. 3(b) and (d) show that 
the primary contributions come from the electrons with a 
broad energy distribution centered at Ep = 2.92 a.u., as 
indicated by a blue arrow. This energy value also agrees 
with the energy relation uj = Ep + Ip, but is far from 
the center of the dip in the red curve in Fig. 2(c). Then 
we can expect that two-center interference has a trivial 
influence on this order, and the first excited state plays 
a negligible role. In Fig. 3(b) of the molecular case, 
the black curve almost superposes the red curve and this 
order in Fig. 1(a) is not suppressed. 

From Fig. 3(b)-(d), we see that around the blue ar- 
rows, there exists a wide energy regime where the in- 
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tegrated amplitudes increase rapidly. For example, for 
the 69th order, the regime spreads from Ep = 2.5 a.u 
to 3.5 a.u., which converts to the energy of about 17 
field-photons. This reveals an important result that con- 
tinuum electrons with a broad energy distribution con- 
tribute equally to one harmonic. The underlying physics 
is that, from the quantum viewpoint, the recombina- 
tion electrons probably emit and absorb additional en- 
ergy from the laser field in the vicinity of the cores p^. 
In particular, the blue arrows here indicate the energy 
Ep = Lu ~ Ip. This energy corresponds to the semiclas- 
sical energy prediction of i?[k(i)] = uj — Ip, obtained 
in the SFA for the instant momentum k(i)[^]. There- 
fore, our simulations presented here provide a quantita- 
tive comparison between the full-quantum energy distri- 
bution and the SFA. The comparison gives significant 
suggestions for proper understandings of the HHG in a 
full-quantum picture, especially for the intramolecular 
interference effect, as shown in Fig. 3(a). Our further 
analyses show that the energy distribution is broader for 
lower harmonic orders, where the Coulomb potential has 
a more important role. We expect that the Coulomb ef- 
fects are mainly responsible for the physics behind our 
simulations. The differences between the SFA and the 
TDSE have been reported in some references [H Him. 
In addition, since the center of the broad energy distri- 
bution locates at Ep = u> — Ip, it confirms the relation 
LD = Ep + /p, which holds in an averaging process as 
shown in Ref.[i,[ll. 

From the above discussions, we anticipate that the two- 
center interference effect in the molecular HHC not only 
relies on the population of the first excited state, but also 
could be affected by laser parameters such as intensity. 
This is shown in Fig. 1(c). First, the plottings in Fig. 
1(c) verify the important role of the first excited state 
in the molecular HHG. Furthermore, the comparison be- 
tween the black curves in Fig. 1(a) and (c) reveals that 
the pronounced minimum in the striking hollow region 
shifts as the laser intensity changes. Specifically, it shifts 
from the 27th order in Fig. 1(a) to the 23rd order in Fig. 
1(c), as indicated by the vertical arrows. It should be 
mentioned that in Ref.0], a spectrum-smoothing proce- 
dure is used to help the identification of the interference 
minima. This procedure is not adopted in our analysis, 
since we expect a "direct" comparison of the numerical 
observations and the experimental measurements. 

The above interference picture is somewhat different 
from the point-emitter picture. If we consider that only 
the electrons, that transit back to the ground state and 
have the energy Ep — u) — Ip, contribute significantly 
to one harmonic wQ, Eq. 2 can be approximated to 
F{lo) = (0|Vy|p)/dtaS(t)e*"*Cp(t), with to = Ep + Ip. 
The above expression explicitly shows that the inter- 
ference minimum in the molecular HHG spectrum that 
corresponds to the minimal extremum of the transition 
dipole (0|Vy|p) is independent of the laser intensity. 

In the full expression of Eq. 2, it says that the 
electrons, that transit back to varied bound states and 



have diverse energy, contribute to one harmonic. Our 
ah initio calculations reveal that the primary contribu- 
tions to one harmonic come from the recollision elec- 
trons with a broad energy distribution, each of which 
is weighted by the probability amplitude G(p, u) = 
E„'("'|VF|p) / dtal,{t)e''^^Cp{t). This amplitude de- 
pends on the laser intensity through the terms a„' {t) and 
Cp{t). The first excited state also has a nonnegligible role 
in the molecular HHG. Therefore, the interference mini- 
mum in the molecular HHG spectrum could be regulated 
by the laser intensity. We noted that the cutoff in Fig. 2 
of Ref. 3 is probably at the 45th order that corresponds 
to the laser intensity of / = 3 x 10^"* W/cm^. Conipared 
to the intensity of / = 2x lO^'' W/cm^ used in Ref.Q, we 
expect that the uncertain in the calibration of the laser 
intensity can occur there and the different observations 
in Ref. 0, Q can be attributed to the dependence of the 
interference pattern on the laser intensity. 

We now extend our considerations to the 2D and 3D 
cases. For the 2D H^, diagonalizing Hq becomes compu- 
tationally intractable, and instead of Eq. 2, we simulate 
using F{uj) ^ /E„' <'(0("'|e • Vy|^/'(i))e^'^*dt. For 
the 3D hydrogen-like atom, we simulate using F{uj) — 

/dpG'(p,w) with G'(p,w) = Y.Z'~f '^n',Lm{V-,^) and 
G„',z,™(p,c^) = (n',;,m|e- VF|p) j'dtal, i ,^{t)e'^'cp{t). 
Here, |p) is the accurate continuum state with outgoing 
wave boundary conditions[2l|,[22|. The definition of n' is 
the same as in the ID case. These results are presented in 
Fig. l(d)-(f). For the 3D atom. Fig. 1(d) shows the ap- 
plicability of our numerical scheme. Analysis on the inte- 
grated contribution of | X]n'7° /o' ^P^"'.'.™(P' 
obtains the similar results as in the ID case. We plot 
the results in the right panels of Fig. 3. For the 2D 
molecule. Fig. 1(e) and (f) confirm the important role of 
the excited state in the hollow region (as indicated by the 
dotted lines) and the shift of the pronounced minimum 
(from the 39th to the 51st orders as indicated by the solid 
lines) as the laser intensity changes. These major results 
are also validated by our further simulations for H^ with 
varied internuclear distances. Note that this shift arises 
from not only the nontrivial role of the excited state in 
the molecular HHG but also the fact that continuum elec- 
trons in a broad energy regime contribute importantly to 
one harmonic. Even we only consider the transition back 
to the ground state, the corresponding green curves in 
the molecular cases in Fig. 1 also reveal the shift of the 
pronounced minimum. Particularly, all of our ID and 2D 
results show that the first excited state already leaves its 
unambiguous footprints in the molecular HHG spectrum. 
This finding throws a new light in the difference of the 
spectral amplitudes for the molecule and its reference 
atom in the orbital tomography procedure 

In summary, we have studied the HHG with a full- 
quantum treatment. Our simulations revealed a wide 
regime of continuum electronic energy that contributes 
significantly to one harmonic, and a nonnegligible con- 
tribution of the first excited state to the molecular HHG. 
Consequently, the two-center interference pattern is in- 
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flucnccd by the laser intensity, as consistent with exper- 
iments. These results have important implications for 
ultrafast imaging of transient molecules. 
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